Unified Scaling Law for Earthquakes 
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We show that the distribution of waiting times between earthquakes occurring in California obeys 
a simple unified scaling law valid from tens of seconds to tens of years, see Eq. (1) and Fig. 
4. The short time clustering, commonly referred to as aftershocks, is nothing but the short time 
limit of the general hierarchical properties of earthquakes. There is no unique operational way of 
distinguishing between main shocks and aftershocks. In the unified law, the Gutenberg-Richter b- 
value, the exponent —1 of the Omori law for aftershocks, and the fractal dimension df of earthquakes 
appear as critical indices. 
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Earthquakes are a complicated spatio-temporal phe- 
nomenon. The number of earthquakes with a magnitude 
M > to is given by the Gutenberg-Richter law |jj . In ad- 
dition to the regularity in the rate of occurrence, earth- 
quakes display a complex spatio-temporal behavior . 
The spatial distribution of epicentres is fractal and they 
occur on a fractal- like structure of faults Short- 
range temporal correlations between earthquakes are ex- 
pressed by Omori's law Jj|, which states that immediately 
following a main earthquake there is a sequence of after- 
shocks whose frequency decays with time as T~ a , a ~ 1. 
This has led to the commonly held belief that aftershocks 
are caused by a different relaxation mechanism than the 
main shocks. 

The observed temporal complex behavior is obviously 
of dynamical origin. However, the statistics of earth- 
quakes as well as the geometrical fractal structure dis- 
played by the faults and by the spatial distribution of 
epicenters is also a result of a dynamical process and one 
might speculate whether it is possible to unify these ob- 
servations. 

We propose a unified scaling law for the waiting times 
between earthquakes, expressing a hierarchical organiza- 
tion in time, space, and magnitude. There is a correlated 
regime where the distribution of waiting times between 
earthquakes is a power-law, T~ a ,a rj 1 and an uncor- 
rected regime. However, the waiting time interval for 
the crossover between the two regimes for earthquakes 
larger than a given magnitude depends on the area and 
magnitude under consideration. 

An earthquake catalogue covering the period 1984 — 
2000 in a region of California spanning 20°N - 45°N lat- 
itude and 100°W - 125°W longitude was analyzed §. 
The total number of recorded earthquakes in the cata- 
logue is 335, 076. The number of earthquakes N(M > to) 
with magnitude larger than to is given by the Gutenberg- 
Richter law H log w N(M > m) oc -bm, b w 0.95 see Fig. 
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Magnitude m = log 1(l (S) 

FIG. 1. The number of earthquakes N(M > m) 
with a magnitude larger than m per year (open cir- 
cles). The dashed line is the Gutenberg-Richter law 
log 10 N(M > m) oc —bm, b — 0.95. The deficit at small mag- 
nitude m < 2, is related to the problems with detecting small 
earthquakes, so only earthquakes with m > 2 will be consid- 
ered. 

The spatio-temporal analysis was carried out as fol- 
lows. We covered the region with a grid with cells of 
size L x L, see Fig. |2|and defined the waiting time T as 
the time interval between the beginning of two succes- 
sive earthquakes. We then measured Ps,l(T), the distri- 
bution of waiting times T, between earthquakes occur- 
ring within range L whose magnitudes are greater than 
to = log (5) . 
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FIG. 2. The region of California divided into grids of two different cell sizes (a) L = 4° and (b) L = 1°. 



Figure || shows the resulting set of curves Ps,l(T), for 
time scales ranging from seconds to 16 years, for several 
values of S and L, plotted on double logarithmic scale. 
Obviously, the curves differ widely. Some general trends 
can be seen, however. There is a linear regime, indicat- 
ing a power-law distribution, extending up to a cut-off 
indicating an upper limit of the waiting time. For fixed 
cell size L and increasing cut-off S (or m), the range of 
the power-law regime increases. For fixed cut-off S and 
increasing cell size L, the range of the power-law regime 
decreases 
l(T 




FIG. 3. The distribution Ps,l(T) of inter-occurrence 
time T between earthquakes with magnitude m greater 
than logio(S) within an area of linear size L. The 
solid circles, squares, and triangles correspond to cut-offs 
m = log 10 (5') = 2,3 and 4, respectively. The color coding 
represents the linear size L = 0.25° (black), 0.5° (red), 1° 
(green), 2° (blue), and 4° (orange) of the cells covering Cal- 
ifornia. For T < 40 s, earthquakes overlap and individual 
earthquakes cannot be resolved. This causes the deficit for 
small T, so only intervals T > 38 s will be considered in 
the following. Notice, that for fixed cut-off magnitude m but 
decreasing linear size L, the deviation from the Omori law 
T~ a ,a m 1 sets in for larger values of T. On the contrary, 
with fixed linear size L but decreasing cut-off magnitude m, 
the deviation from the Omori law T~ a ,a ps 1 sets in for 
smaller values of T. 

In Fig. (§ the curves are re-plotted in terms of re- 
scaled coordinates. The x-axis is chosen as x — TS~ b L d f , 
and the y-axis represents y = T a Ps.L(T). The rescaling 
causes a shift of the curves in Fig. || that depends on L 
and S. For suitable choice of the interval exponent a, the 
magnitude exponent b, and the spatial dimension df, all 
the data collapse nicely onto a single well-defined curve 
f(x), that is, 



T a P s , L (T) = f{TS-"L d f) 
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FIG. 4. The data in Fig. g with T > 38 s 
re-plotted with T a Ps,l(T) as a function of the variable 
x = cTS~ b L df , c = 10~ 4 . The data-collapse implies a unified 
law for earthquakes. The Omori law exponent a = 1, Guten- 
berg-Richter value b = 1, and fractal dimension df — 1.2 have 
been used in order to collapse all the data onto a single, unique 
curve f(x). The estimated uncertainty in the exponents is less 
than 0.2. The function / is constant for x < 1, corresponding 
to the correlated Omori law regime, while it is decaying fast 
for large arguments x > 1, associated with the uncorrelated 
regime of earthquakes. Whether two earthquakes are to be 
categorized as belonging to a correlated or uncorrelated se- 
quence not depend independently on the values of T, S and L 
but only on the value of the product x = TS~ b L df . 

This equation expresses the unified scaling law for 
earthquakes. The function f{x) consists of a constant 
part and a decaying part, separated by a sharp kink. 
The constant part corresponds to the linear, power-law 
part in Fig. [| since we have multiplied Ps,l(T) with 
T a . Any deviation from power-law behavior would show 
up dramatically in this type of plot. Nevertheless, the 
function is approximately constant over 8 orders of mag- 
nitude! The rapidly decaying part is consistent with an 
exponential decaying function implying an uncorrelated 
regime for large values of x. This is indeed what one 
would expect on physical ground: earthquakes that are 
separated by large enough distances or long waiting times 
will be uncorrelated. 

The index a « 1 can be identified as the Omori-law 
exponent for aftershocks, b « 1 is the b- value in the 
Gutenberg-Richter law, and df « 1.2 describes the 2d 
fractal dimension of the location of epi-centers projected 
onto the surface of the Earth. 

The data collapse implies that the waiting time dis- 
tribution depends on T, S, and L only through the vari- 
able x. Only critical processes exhibit this type of data 
collapse, known as scaling in critical phenomena [||, so 
our analysis demonstrates that earthquakes are a self- 
organized critical (SOC) phenomenon p-pTj], as had been 
anticipated from the existence of the Gutenberg-Richter 
law p^-p^[. The data collapse shows that there is no 
separate relaxation mechanism for aftershocks. The three 
exponents a, 6, and df characterizing earthquakes emerge 



as critical indices in the unified law. The estimated un- 
certainty in the values of the critical indices is less than 
0.2. Whether the critical exponents vary with region and 
maybe even with time is an interesting question that is 
outside the scope of the present letter but we urge further 
studies in that direction. 

Depending on the value of scaling argument x, and 
thus the chosen values of L and m (or S) , two successive 
earthquakes will either be correlated, for x small, (i.e. to 
the left of the kink in Fig. |]), or uncorrelated, for x large 
(i.e. to the right of the kink in Fig. ^). 

Depending on the length scale L of observations, and 
the magnitude magnitude m (or S) chosen, the corre- 
lated Omori T~ a regime may range from seconds to tens 
of years (and probably much longer if data were avail- 
able). If the earthquakes are correlated they may be in- 
terpreted as belonging to an aftershock sequence. If they 
are uncorrelated, they may be interpreted as main events. 
However, this depends on L and m through the variable 
x = TS~ b L df and has no absolute meaning. Therefore, 
there is no unique way of characterizing earthquakes as 
aftershocks or main events. 

To summarize, the short time correlations given by 
Omori's law is just the short time limit of a general hi- 
erarchical scaling phenomenon occurring at all accessible 
time scales. Amazingly, the statistics of aftershocks oc- 
curring within minutes of an earthquake can be simply 
related to the statistics of earthquakes separated by tens 
of years! 

One may think of the value of m (or S) at the kink as 
being a 'characteristic' magnitude of earthquakes. How- 
ever, since the position of the kink is a function of x, not 
to (or S), one cannot specify the magnitude of the 'char- 
acteristic' earthquake without at the same time specify- 
ing a time scale and an area. However, there is no special 
time scale that can play any absolute role for the dynam- 
ics of earthquakes, limited at the upper end by the time 
scale of tectonic plate motion, and at the lower end by 
the duration of earthquakes. 

How should one physically understand the fundamen- 
tal law, Eq. (1)? Let us first discuss the meaning of the 
scaling variable x = TS~ b L d ' . The quantity S~ b L d f ap- 
pearing in it's definition is a measure of the average num- 
ber of earthquakes per time unit with magnitude greater 
than m = logio(S) occurring within the range L. Thus, 
a; is a measure of the average number of such earthquakes 
occurring within a time interval T. The law states that 
the distribution of waiting times depends only on this 
number. When this number exceeds a well-defined value 
(the position of the kink in Fig. |i|), the earthquakes 
sharply become less correlated. 

Think of earthquakes being generated by 'processes', 
each producing a sequence of correlated earthquakes with 
a T~ a distribution. These processes correspond to a 
sequence of avalanches in self-organized critical models 
of complex phenomena. Visually, one might think of 
the processes as the activity associated with dynami- 
cally changing fault segment patterns. The law indi- 
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cates that the crust operates in the true SOC slow-driving 
regime [ [l6| where the individual processes (avalanches) 
do not overlap. Because of the non-zero driving rate, 
several spatially separated processes are active simultane- 
ously. The kink on the /-curve indicates the point where 
one crosses into the regime where spatially independent 
earthquakes, belonging to different processes, are sam- 
pled within a window spanned by T and L. For small 
enough L and T only a single, correlated process is sam- 
pled. 
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